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Abstract 

The ultrahigh-energy (UHE) proton and neutrino spectra resulting from collapse or an- 
nihilations of topological defects surviving from the GUT era are calculated. Irrespective 
of the specific process under consideration (which determines the overall normalization of 
the spectrum), the UHE proton spectrum always ‘recovers’ at ~ 1.8 x 10 ll GeV after a 
partial Greisen-Zatsepin-Kuz’min ‘cutoff’ at ~ 5 x 10 10 GeV and continues to a GUT-scale 
energy with a universal shape determined by the physics of hadronic jet fragmentation. 
Implications of our results axe discussed. 
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Topological defects[l], e.g., monopoles, cosmic strings, domain walls, superconducting 
cosmic strings, etc., and various possible hybrid objects made of these, are likely to have 

i 

been formed in symmetry-breaking phase transitions in the early universe. Because of 
their topological stability the defects can survive indefinitely, until and unless physically 
destroyed due to collapse, annihilations or other reasons[2-5]. When topological defects are 
destroyed the energy trapped in them is released in the form of massive quanta (hereafter 
generically referred to as X particles) of the various fields (gauge bosons, higgs bosons, 
superheavy fermions) that ‘constitute’ the defects. The X-particles can then decay into 
known quarks, gluons, leptons, etc, which ultimately materialize into, among other par- 
ticles, nucleons, gamma rays and neutrinos with energies up to ~ m x , the mass of the 
X-particles released from the defects. If the defects were originally formed in a phase tran- 
sition at the GUT-scale ~ 10 16 (reV, then we have here a possible natural mechanism of 
production of ultrahigh-energy (UHE) cosmic ray (CR) particles up to an energy of this 
order, without any acceleration mechanism. Below we give a general calculation of the 
expected evolved proton as well as neutrino spectra resulting from this kind of processes. 
We assume a = 1 “flat” universe and a Hubble constant H 0 = 100. h Km.s^ 1 .Mpc~ l , 
with h = 0.75 throughout. 

The rate of release of X-particles due to destruction of topological defects can, in 
general, be effectively expressed in terms of two fundamental parameters entering in the 
problem, namely, the mass scale m x and the Hubble time ~ t, in the form (with % — c = 


Mpiy/G = kg = 1 ) 


dn x (tj) 

dti 


= 


p +-4+P 


( 1 ) 


where k and p are dimensionless constants whose values depend on the specific process 
involving specific kind of topological defect, (the ‘amplitude’ « may, in general, depend on 
p ), and n x {ti ) denotes the number density of the X-particles released at the time t{. For 
example, p = 1 for a system of collapsing cosmic string loops[4,5] as well as for a system of 
collapsing monopolonia[2,7], while p = 0 for a system of saturated superconducting cosmic 
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string loops[3]. We shall see below that the observed UHE CR flux gives an upper limit 
to the possible value of k for any given value of p. 

We assume that each X decays into a quark and a lepton each carrying an energy m x / 2. 
The quarks fragment and produce jets of hadrons. We use the following hadronic frag- 
mentation formula[2,3] that is consistent with leading log QCD behavior and reproduces 
well the particle multiplicity seen in GeV-TeV jets in colliders: 


~ 0.08 exp (2.6y/ln(l/x)j (1 — x) 2 (x'y/ln( l/r)j . ( 2 ) 

In eq.(2), x = E being the energy of any hadron in the jet, and Ej et ^ m x / 2 

is the total energy in the jet, and N h is the number of hadrons carrying a fraction x of 
the total energy in a jet. It turns out[5] that, in the range 10~ 10 < x < 10 -2 , Eq.2 is 
numerically well approximated by a power-law of the form 


dN h 

dx 


~ 3.322 r -1324 , 10" 10 < x < 10~ 2 , 


( 3 ) 


which gives rise to a considerably flatter injection spectrum (see below) than the usually 
assumed UHE CR injection spectrum having a power-law index ~ —2.5. 

Most of the hadrons in a jet will be pions, but a small fraction (~ 3%) will be nucleons 
and antinucleons. For the propagation times and energies of our interest most (but not all) 
of the neutrons and antineutrons in a jet will /?-decay to protons and antiprotons which 
we shall hereafter collectively refer to as protons. From Eq.l, the injection spectrum, 
$>p(Ei, ti ), for protons, i.e., the number density of protons injected per unit energy interval 
at an injection energy Ei per unit time at an injection time is (with x = 2 Ei/m x ) 

X 0.03 X 

dti m x dx 

Each charged pion (in the jets) after decay ultimately gives rise to three neutrinos — one 
v pi one z/ M , and one u e (or P e ) — all of roughly the same energy. The z/ M + injection 
spectrum resulting from charged pion decay can be written as[8] 

$JEi,ti) = 2 x ^i(ti) x 2.343 x — x 0.97 / ~^±dx. 

dt i m X J 4 .686Ei/m x x dx 


x dx 



Note that, except for the overall source evolution factor dn x /dti given by Eq. 1, the basic 
form of injection spectrum, (for both protons and neutrinos) is determined by the micro- 
physics of hadronic jet fragmentation and not by any external astrophysical parameters. 

For any injection spectrum $(Ej,t,) due to a isotropic source distribution, the particle 
flux (i.e., the number of particles crossing unit area per . unit time per unit solid angle per 
unit energy interval) today, for ultrarelativistic particles traveling along straight-line paths 
in a — 1 universe is 




E 0 


3 f ZiT 


.(J5o) 


dzif(Ei t £,•), 


where 




^ „ / 1 ,7372xl0 19 aec 


U(- 


E 0 


if z , < z tq , 
if z ,• > z eq . 


■) gU 3 s: 1/2 (i + *0 (M) 

In the above equations z,- is the redshift corresponding to the injection time t t , z eq is 
the redshift of equal matter and radiation energy densities, t 0 is the present age of the 
universe, Eq is the energy of the particles in the present epoch, E t - = Ei(Eo,Zi) is the 
injection energy, Zi tmax (Eo ) is the maximum possible injection redshift corresponding to 
the earliest possible injection time ti tm i n {E o) for a given Eo, and and g m are the total 
number of effective relativistic degrees of freedom contributing respectively to the entropy 
density and energy density of the universe[9] at the epoch z{. The upper limit on zi jTnax is 
given by Ei(E 0 ,z ijmax (E 0 )) = m x j 2. 

The energy-loss processes for UHE protons propagating through the cosmic medium, 
including the important Greisen-Zatsepin-Kuz’min (GZK)[lO]-effect — whereby UHE pro- 
tons above a certain threshold energy (~ 5 x 10 10 GeV) catastrophically lose energy due 
to photopion production off the cosmic background photons — have been studied by many 
authors; for comprehensive discussions see, for example, Refs. [11,12]. For UHE protons 
we obtain the function E{(Eq, z,) by numerically solving the differential equation for energy 
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loss discussed in Ref. [12]. We find that for UHE protons Z{ %max < 1.25 at all energies of 
our interest. 

For UHE neutrinos the relevant process is the absorption[13] of neutrinos through the 
process v -f where i — e, p, r, it, d, s, c, etc., and i>u represents the thermal 

background neutrinos (present temperature 1.9 °K). Denoting by z a (E Uj o) the redshift 
above which neutrinos of present energy E „ t o are absorbed due to the process mentioned 
above, we get 


1 + z a (E u ,o) 


/ _ 1/3 

8.8143 x 10 s , for E„, 0 < 3.08 x 1 0 5 GeV, 

4.8286 x 10 5 , for E„ l0 > 3.08 x 10 5 GeK 


Finally, for neutrinos, we have 


, v . / . 0.2134 m x 

1 4- Zi,max(E»,o ) = mm (1 + z a (E „ t 0 ), — „ 

where we have used the fact[14] that the maximum neutrino energy is 0.4268 times the 
maximum pion energy (~ m x /2). 

The proton and neutrino differential spectra obtained by taking into account the energy 
loss of protons and absorption of neutrinos are shown in Figures 1 and 2 respectively for 
m x = 10 15 GeV. For any given value of p (see Eq. 1) we obtain an upper limit (« m ai) on 
the ‘amplitude’ k for the process by requiring that the calculated UHE proton flux at any 
energy not exceed the observed UHE CR flux[15]; the results shown axe for « = K max . The 
evolved proton spectrum due to a proton injection spectrum oc (1 + Z{) 9 / 2 E~ 2 5 (hereafter 
called the “test spectrum”), continued to the energy m x / 2, is also shown in Fig. 1 for 
comparison. (Note that the index 9/2 of the redshift factor corresponds to the case of p = 1 
in Eq. 1). Comparison of the short-dashed curves (the ‘redshifted’ injection spectra) with 
the solid curves in Fig. 1 shows the effect of ‘photopion energy loss’ (the GZK-effect) on the 
UHE proton flux — the so-called “pile up” or “bump”[ll,12] at Eq ~ 5 x 10 10 GeV is clearly 
visible. Contrary to usual assumption, however, there is really no complete GZK-‘cutoff’ 
of the UHE proton spectrum associated with the photopion energy loss; instead Elj(Eo) 
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only takes a “dip” after the “pile up” and then rises again as Eq increases further until 
it turns over before the final cutoff at E 0 ~ m x /2. The spectrum flattening (‘recovery’) 
at Eq ~ 1.8 X \O xl GeV after the GZK-‘cutoff’ at Eq ~ 5 x 10 10 GeV is partly due to a 
progressive decrease (with increasing energy) of the slope of the ‘photopion energy loss’ 
curve[l2] for protons beginning at Eq ~ lQ lx GeV, and, of course, due to continuation of the 
injection spectrum to a much higher energy. Note that the recovery of the test spectrum 
(which, properly normalized, reproduces the observed UHE CR spectrum reasonably well 
at the observed energies; see Fig. 1), is much weaker because of its steeper injection 
spectrum. 

The protons surviving at energies above the GZK-‘cutoff’ originate at small injection 
redshifts, and so the effect of source evolution (specified by the value of p) on the final 
shape of the proton spectrum at these energies is small compared to that due to photopion 
energy loss. This, and the fact that the shape (though not the overall normalization) of 
the injection spectrum is determined by the physics of hadronic jet fragmentation inde- 
pendently of the value of p, lead to a universal (independent of the value of p) spectrum 
shape for protons at these energies. The shape of the UHE neutrino spectrum, on the 
other hand, reflects the effect of source evolution (i.e., dependence on p ) because of higher 
possible injection redshifts. 

For neutrinos we see from Fig. 2 that the lowest energy above which neutrinos from 
topological defects may be ‘visible’ above the atmospheric neutrino background[l6] is 
smaller for processes with smaller values of p. 

From Fig.l it seems unlikely that topological-defect-induced UHE particles give any sig- 
nificant contribution to the observed UHE CRfiux at energies below ~ 5xl0 10 GeV r . Above 
this energy, however, the topological defects can in principle give significant contribution 
to any (observed?[l7]) CR flux. So, detection of UHE CR events at Eq » 5 x 10 10 GeV 
may give us a signature of GUT. in general and topological defects in particular. From 
Fig. 1 a rough power-law approximation to our universal UHE proton spectrum above 
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the GZK-‘cutoff’ gives j(Eo) oc E 0 3-8 for 4.5 x lO^GeV < Eq < 1.8 x 10 n GeV and 
j(E 0 ) oc Eq 1 * 3 for 1.8 x 10 n GeF < E 0 < 5 x 10 13 GeF. 

The maximal expected integral flux for both protons and neutrinos axe shown in 
Fig. 3. The maximal proton integral flux above 7 x 10 10 GeV, for example, is ~ 
5 x 10 -16 m -2 .s“ 1 .sr“ 1 giving an event rate of ~ 5 x 10 ~ 2 yr~* .KmT 2 in a 7r steradian 
detector. The numbers for neutrinos depend on the value of p and can be read off from 
Fig. 3. Thus for processes with values of « not too small compared to the maximal values 
obtained above, the resulting UHE proton flux may be detectable in some of the planned 
large Extended Air Shower Arrays such as[18] the 100 Km 2 AGASA experiment in Akeno, 
Japan, the “HiRes” detectors of the Fly’s Eye experiment, and the proposed lOOOifm 2 
array in USSR. Similarly the UHE neutrinos from these processes may be detectable in 
experiments like DUMAND[19]. 

To summarize, characteristic features of topological defect-induced hadronic CR spec- 
trum include: (1) A spectrum flattening (‘recovery’ ) at ~ 1.8 x 10 11 GeU after the GZK 
steepening (‘cutoff’ ) at ~ 5 x l0 lo GeU, and continuation of the spectrum to a GUT scale 
energy with a universal shape ; (2) hadronic component consists of protons and perhaps 
some neutrons, but no “heavies” such as Fe nuclei; and (3) presence of antiparticles in 
the primary spectrum, which may however be difficult to distinguish from particles at the 
extremely high energies we are considering. 

Finally, the numerics! results given above are for a specific value of m x = 10 15 GeV. 
The flux of particles at any given energy will be higher for a smaller value of m x , but 
then the final cutoff of the spectrum shifts to a smaller value such that overall energy 
conservation is maintained. 

PB is supported by NSF grant INT 91-10117 to the University of Chicago, CTH by 
DoE at Fermilab, and DNS in part by NSF grant AST 91-22629, by DoE grant DE-FG02- 
91 ER 40606 and by NASA grant NAGW 1321 at University of Chicago and by the DoE 
and NASA grant NAGW 1340 at Fermilab. 
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Figure Captions: 

Fig. 1. The UHE proton spectrum due to topological defects. The four solid curves 
(which are indistinguishable for Eq > 5 x 10 10 GeV) are the maximal flux spectra for four 
different values of p. These correspond, from the lower-most solid curve to the uppermost 
solid curve, to p = 1.5, 1, 0.5 and 0, respectively, the corresponding values of K max being 
1.18 x 10 -29 , 0.2313, 4.52 x 10 27 and 8.82 x 10 55 , respectively. The four short-dashed 
curves (“redshifted injection spectra”) correspond to the same four values and order of 
p and corresponding K max as above, but now calculated with no other energy loss of 
protons except that due to expansion of the universe. The dashed curve (terminating at 
7 x 10 1Q GeV) represents the UHE CR flux given by the Fly’s Eye group[15], and the dash- 
dotted curve immediately above it is for the “test spectrum” (see text). The normalization 
of the last curve has been conveniently chosen for clarity of display. 

Fig. 2. The UHE neutrino spectra. The solid curves represent the maximal flux 

for the same set of values and order of p and corresponding K. max as in Fig. 1.. The dashed 
curve represents the estimated horizontal atmospheric neutrino fiux[l6]. 

Fig. 3. The maximal integral UHE proton (p) and neutrino (u) spectra for the same 
values and order of p and corresponding « mox as in Figs. 1 and 2. The dashed curve 
represents the UHE CR flux given by Fly’s Eye group[l5] and the dash-dotted curve 
immediately above it is for the “test spectrum” with the same normalization as in Fig. 1. 
The two short-dashed curves marked vD and i/U represent the upper limits on downward- 
going and upward- going neutrinos, respectively, given by the Fly’s Eye Group[20]. 
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